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r- 1 . Abstract 

We investigate possible search strategies for QCD-instanton induced processes at HERA 
in the deep-inelastic scattering (DIS) regime. Our study is based on the Monte Carlo 
_C ■ generator QCDINS for instanton-induced events and the standard generators for normal DIS 

events. It appears possible to isolate an instanton enriched data sample via an optimized 
multi-dimensional cut scenario for a set of six most instanton-sensitive DIS observables. 
As a further central point, we investigate the stability of our results with respect to a 
variation of the (hadronization) models available for the simulation of both normal DIS and 
instanton-induced events. Dependencies on the variation of certain inherent parameters are 
also studied. Within the "bandwidth" of variations considered, we find that the normal DIS 
background is typically much more sensitive to model variations than the I-induced signal. 



* Contribution to the Proceedings of the DESY Workshop 1998/99 on Monte Carlo Generators for HERA 
Physics. 



1 Introduction 

Instantons (I) have been known to theoretical physics since the mid seventies [|l| |2|, || ||. They 
represent non-perturbative fluctuations of the gauge fields in non-Abelian theories like QCD and 
weak interactions, associated with tunneling transitions between degenerate ground states (vacua) 
of different topology 

Instanton transitions induce processes |2], |3|, ||] violating the (classical) conservation of certain 
fermionic quantum numbers, namely chirality in QCD and the sum of baryon and lepton number 
in weak interactions. These processes are forbidden in ordinary perturbation theory, but have to 
exist in general due to the famous AB J chiral anomaly J5], ^|, [7] . 

The DIS regime at HERA offers a unique possible discovery window for events induced by QCD- 
instantons through their characteristic final state signature || |S], [II], O, [T2j and a sizable rate, 
calculable^ within I-perturbation theory [14], [15|, [16j . An experimental detection of these processes 



would correspond to a novel, non-perturbative manifestation of non-Abelian gauge theories and 
would clearly be of basic significance. 

This report is organized as follows^: 

We start off in section |^ by introducing the dominant I-induced DIS process along with the 
relevant kinematics. Then we proceed by summarizing the essential ingredients and the basic 
structure of the Monte Carlo generator QCDINS for instanton-induced DIS events || |18| , on 
which the present study is based. In section [|, we recall the theoretical results for the I-induced 
cross section at HERA and summarize the corresponding characteristic topology of the hadronic 
final state. Section [| is devoted to a study of the possibility to reconstruct the Bjorken variables 
of the I-subprocess. Sections | and ^| contain the central results of this investigation. In Section |5], 
we report on a possible search strategy for these processes. The goal is to isolate an instanton 
enriched data sample via an optimized multi-dimensional cut scenario for a set of most I-sensitive 
DIS observables. Finally, in section || we investigate the stability of our results with respect to 
a variation of the (hadronization) models available for the simulation of both normal DIS and 
instanton-induced events. Dependencies on the variation of certain inherent parameters are also 
considered. 

All the studies presented in this report are performed in the hadronic center of mass frame?] 
(hCMS), which is a suitable frame of reference in view of a good distinction between I-induced 
and normal DIS events (c. f. [[H]]). The results that we obtain are based on a study of the hadronic 
final state, with typical acceptance cuts of a HERA detector being applied in the laboratory frame 
(15° < Qhadron < 155° and pT(hadron) > 0.15 GeV, for charged particles, Qhadron > 4°, for all 
stable particles, and 0.25 < p^(K°) < 4.5 GeV 2 , for neutral kaons)f]. 

^Tor an exploratory calculation of the I-contribution to the gluon structure function see Ref. jig] . 
2 An extended version of this paper will appear elsewhere []l7j . 

3 This frame is defined by q + P = 0. Therefore, the photon and the incoming proton collide head-on. 
4 In the laboratory frame the incoming proton points in the positive z-direction, while in the hCMS the proton 
points in the negative z-direction. 



2 The Monte Carlo Generator QCDINS 



In deep-inelastic e ± P scattering, I-induced events are predominantly associated with a process 
structure as sketched in Fig. [J]: A photon splitting into a quark- ant iquark pair, fuses with a gluon 
from the proton in the background of an instanton or an anti-instanton. Besides the current- 
quark (jet), the partonic final state consists of 2n/ — 1 right (left)-handed (massless) quarks 
and anti-quarks and n g gluons emitted by the instanton (anti-instanton). Correspondingly, the 
total violation of chirality = ^Qr ~ #?l is ±2n/, in agreement with the ABJ chiral anomaly 
relation @, |, 0. 

DIS variables: 

S = (e + P) 2 

Q 2 = -q 2 = -(e - el) 2 

x Bj = Q 2 I (2P ■ q) 

Z/Bj = Q 2 I (Sx B i) 
W 2 = ( q + pf = Q2 ( 

s = {q + gf 

i = x Bi (1 + s/Q 2 ) 
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Figure 1: Kinematic variables of the dominant I-induced process in deep-inelastic scattering. 
I-induced processes initiated by a quark from the proton are suppressed by a factor of a 2 with 



respect to the gluon initiated process [|H|. This fact together with the high gluon density in the 



relevant kinematical domain at HERA, justifies to neglect quark initiated processes. 



QCDINS H, [18] is a Monte Carlo package for simulating QCD-instanton induced scattering pro- 
cesses in DIS. It is designed as an "add-on" hard process generator interfaced by default to the 
Monte Carlo generator HERWIG PH . 

QCDINS incorporates the essential characteristics, that have been derived theoretically for the 
hadronic final state of I-induced processes: notably, the isotropic production of the partonic final 
state in the I-rest system (q'g CMS in Fig. [I]), flavour "democracy", energy weight factors different 
for gluons and quarks, and a high average multiplicity 2rif + 0(l/a s ) of produced partons with 
a (approximate) Poisson distribution of the gluon multiplicity. 

Let us briefly summarize the main stages involved in QCDINS, to generate the complete I-induced 
partonic final state in DIS. 

The first stage is the generation of the various Bjorken variables (c. f. Fig. |l|) of the I-induced 
DIS process. They are generated in the order Q' 2 , x', £, xbj, Ubj, with weights corresponding to the 



I-induced total cross section [|15|, |16|| , 
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subject to appropriate kinematical restrictions (e.g. < xbj < ^bj/£ < x' < 1) and fiducial cuts. 
Here e 2 , denotes the electric charge squared of the virtual quark q' in units of the electric charge 

squared, e 2 = Ana, and f g is the gluon density in the proton. The weight factor P ' , accounts for 



the flux of virtual quarks q' in the I-background [|15 , |T 
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In Eq. ([!]), the I-induced g'g-subprocess total cross section a , contains the essential instanton 
dynamics ]n| . As illustrated in Fig. |], it is very steeply growing for decreasing values of Q' 2 and 
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Figure 2: I-subprocess cross section fl5|| displayed versus the Bjorken variable Q' 2 with x' fixed (left) 
and versus x' with Q' 2 fixed (right) for rif = 3. The dotted lines indicate the corresponding 
effective /-sizes p* [fm] (left) and //-distances R* in units of p* (right), respectively. 



x', respectively. In order to remain within the region of validity of I-perturbation theory |15, ^0 
the following cuts are implemented by default in QCDINS, 



Q' 2 > Q'Ln = 64 GeV 2 ; 
An additional cut on the photon virtuality, 

SxbjVbj =Q 2 >Q 2 m 



x' > x'i. 



64 GeV 2 



0.35. 



(3) 



(4) 



should in principle be applied in order to warrant sufficient suppression of non-planar contribu 
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tions |14|], which may spoil the validity of Eq. ([!]). The cross section a , depends strongly on the 

4 



QCD scale A^ and weakly on the renormalization scale /i r . By default, these scales are taken 



MS 

as 



A®=0.282GeV; // r = 0.15Q'. (5) 



MS 

Note that strictly speaking, the underlying theoretical framework refers to massless quarks. There- 
fore, the (default) number of flavours was set to n/ = 3 and ALL in Eq. @ was obtained by a 

standard flavour reduction from the DIS average, ALL = 234 MeV in Ref . [^ • The default value 



for /i, r above corresponds to the minimum [jl5|] , da q > g /d[i r = 0. 

In the second stage of the event generation, the 4-momenta g, q, q', q" of the incoming gluon 
g, the virtual photon q, the virtual quark q' and the current quark q", respectively, are filled. 
Sudakov decompositions of these momenta are used to incorporate various constraints, e.g. on 
the momenta squared, g 2 = m 2 {— (0.75 GeV) 2 by default), q 2 = —Q 2 , q' 2 = —Q' 2 , q" 2 = raj?,,. 
The 4-momentum of the outgoing lepton, e', is calculated subsequently. 

In the third stage, the partonic final state of the I-induced g'g-subprocess is generated in the 
q'g CMS as follows. The number n g of produced gluons is generated according to a Poisson 
distribution with mean (n g )^ ~ l/ct s ~ 3 as calculated theoretically in I-perturbation theory. 
n f( = 3) [q ■ ■ -q] - "strings" of partons are set up, each beginning with a quark, followed by a 
random number of gluons and ending with an anti-quark of randomly chosen flavour. There are 
n g + 1 gluons in total and all n/ flavours occur precisely once. A simple example for n g = 3 
outgoing gluons may provide some illustration: [u g g s] [d gu] [s g g d]. A quark and a gluon 
among these 2n/ + n g + 1 partons are (randomly) marked as incoming. 

The momenta Pi of the n = 2nf — 1 + n g outgoing partons in the instanton subprocess CMS are 
then generated according to the energy-weighted phase-space distribution 

(6) 
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corresponding to the leading-order matrix element with different energy weights for gluons and 
quarks. 

Next, the colour and flavour connections of the partons are set up. The colour flow is obtained 
simply by connecting the colour lines of neighbouring partons within each of the above-mentioned 
rif [q . . .q] - "strings" in a planar manner (consistent with the leading order 1/N C expectation). 

The partonic stage of the event generation ends by boosting the momenta of the I-subprocess 
final state partons to the laboratory frame. 

After the perturbative evolution of the generated partons, one may use the hadronization mech- 
anisms inherent in various Monte Carlo models (e. g. HERWIG [19J or JETSET ||22| ) to arrive at 
the complete hadronic final state. 



3 Cross section and Signature 

The total cross section of I-induced events in DIS, calculated within I-perturbation theory, is 
surprisingly large. For xbj > 10~ 3 and 0.1 < y^ < 0.9, the result of Ref. |T3| is 

VjberaW > 0.35; Q'> 8 GeV) ~ 130 pb. (7) 

Hence, given the total integrated luminosity C ~ 30 pb -1 accumulated by each of the HERA 
experiments in the years 1996/1997, in this kinematical region a large number N = c HERA ■ C = 
O(10 4 ) of I-induced events should already have been taken on tape. 

The cross section (|7|) corresponds to a fraction of I-induced to normal DIS (nDIS) events of Jl5| 

a HERA 



However, this prediction for the cross section still contains considerable uncertainties |I5|. One 



QCD scale A^-A in Eq. <^). If Ai± is varied within the allowed range of ±65 MeV, the cross 



of the dominant ones arises from the experimental uncertainty of ±65 MeV associated with the 

fi - Eq. ©. If Ag 
section varies by O(±300)pb and 0(— 100) pb, respectively. 

By far the dominant uncertainty, however, comes in principle from the choice of the lower bounds 
on x' and Q', due to the strong increase of the I-subprocess cross section with decreasing values 
of the cuts x' min and Q' min , respectively (cf. Fig. 0). However, recent high-quality "data" from a 
(non-perturbative) lattice simulation [^] on the topological structure of the QCD vacuum have 
strongly reduced this uncertainty. The lattice results could be directly converted into a "fiducial" 
(x 1 , Q') region, where the predictions from I-perturbation theory should (approximately) hold [20 1 . 



The lower bounds on x' and Q' given in Eq. (0) correspond to this fiducial region and are thus 
expected to be quite reliable . 

Let us next turn to the expected signature of the I-induced final state. 

Fig. |3]a shows the transverse energy distribution in the r)*-<& plane of a "typical" I-induced event in 
DIS, produced by the instanton Monte Carlo generator QCDINS || [Tj|. There are the following 
characteristics: 

• Topology: "Instanton band" and current jet 

In the instanton CMS (the q'g CMS, cf. Fig. |l|) the production of partons coming from 
the I-subprocess is isotropic. This leads to an energy distribution of hadrons restricted to a 
certain range in the pseudo-rapidity r\ (the "instanton-band" ) with a bandwidth of ~ ±1.1 
units in t]. In the instanton CMS this band is localized around rj = 0, while in the hCMS, 
the center of the band is shifted to higher values (depending on the kinematic variables xbj 
and Q 2 ). Besides the band, the hadronic final state of I-induced events exhibits a current 
jet coming from the outgoing current quark (denoted by q" in Fig. III). 
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Figure 3: Transverse energy distributions in the r^-^-plane (hCMS) of the hadronic final state parti- 
cles of an I- induced process (shown for typical detector acceptance cuts), a) "Ideal" event 
signature. Besides the homogeneously populated "instanton-band" (at rf between and 2), 
the current jet at rf ~ 3, 3> ~ 170° is clearly visible, b) The homogeneous distribution of the 
hadrons in the instanton band is destroyed by a current jet with high p?- 

• High multiplicity and "Flavour-democracy" 

In every I-induced event, one pair of quark and anti-quark of all n/(= 3) flavours is simulta- 
neously produced. In addition, a mean of (n g )^ ~ l/a s ~ 3 gluons is expected. Hence, for 
the phase space studied here, we find a mean number of partons of O(10), leading to a mean 
multiplicity of charged particles of O(20) in every event. The actual number of hadrons 
produced mainly depends on the center of mass energy W% accessible. Moreover, due to the 
democratic production of all (effectively massless) flavours, more mesons containing heavy 
quarks should be found in I-induced events than in normal DIS events. The detection of 
K° mesons is experimentally most promising to exploit this feature. 

• High transverse energy 

While normal DIS events exhibit a mean value of 2 GeV of transverse energy per unit of 
|24}1 , for I-induced events a value of the order of 5 GeV in an rf range of < rj* < 4 is 



V 

expected. 



The ideal event topology described above (see also Fig. ||a) is quite often modified by simple 
kinematic effects. In the hCMS the sum of the transverse energy of the particles emerging from 
the I-subprocess and the pt of the current jet has to be approximately balanced. In many cases 



this leads to a destruction of the homogeneous $ distribution of the particles in the band. Fig. |]b 
gives an example for such an event topology. In general, with rising p? of the current jet and 
lower particle multiplicities emerging from the I-subprocess, the band structure becomes less 
clearly visible. 

In addition, the pt of the jet strongly influences the relative if position of the jet and the instanton 
band. Fig. [|a demonstrates this effect for the parton final state. The mean r\* value of the current 
quark increases with rising p? of the current quark, while the mean 77* values of the instanton 
band, outlined by the maximal and the mean {Et weighted) 77* value of the partons coming from 
the I-subprocess, are decreasing. For hadrons (Fig. |]b) the 77* position of the current jet (which 
is reconstructed by a jet-finding algorithm, see section |]) and the Et weighted mean 77* value of 
all particles not belonging to the jet (which we use as an estimator for the central position of the 
instanton band) show qualitatively the same behaviour^. 
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Figure 4: Evolution of the relative 77* position of the partons belonging to the I-subprocess and the 
current jet as a function of the pt of the current quark and the current jet, respectively. In 
a) the instanton-band is outlined by the maximal and the Et weighted mean 77* value of the 
partons coming from the I-subprocess. In b) the position of the band is characterized by the 
Et weighted mean 77* value (calculated without the hadrons belonging to the jet). 



In Fig. |] and throughout this report, the kinematical range given in the context of Eq. (0) has 
actually been extended down to xbj = 10~ 4 with Q 2 = SxbjI/bj > 5 GeV 2 . On the one hand, 
within this enlarged kinematical region, there is the possible influence of non-planar graphs, 
which are not implemented in QCDINS. On the other hand, the hadronic final state topology 
of I-induced events will presumably not change dramatically, being mainly determined by the 
available phase space (see section ||). Since the cross section of I-induced events increases with 



5 The mean 77* values of the centre of the instanton band are lower than those for the parton final state, because 
for hadrons, particles of the proton remnant are also taken into account. Additionally, the mean 77* values for 
PT{current quark) < 3 GeV arc not adequately reproduced on hadron level. In this range, the jet-finding algorithm 
mostly finds the "jet" within the hadrons belonging to the instanton band. 



decreasing values of xbj and Q 2 (as for normal DIS events), it is an experimental challenge to 
investigate the signal to background ratio also in this kinematical domain. In our extended study 
in Ref. |17| we shall restrict ourselves to the "fiducial" high-Q 2 regime (Eq. [|). 



4 Reconstructing the kinematics of the I-subprocess 



In this section we give a short summary of the possibility to reconstruct the kinematic variables 
of the I-subprocess which have been introduced in Fig. |]. A more detailed investigation can 
be found in [P^| . The cross section of the I-subprocess is most sensitive to the values of Q' 
and x'. Therefore the possibility to reconstruct theses variables would allow crucial tests of the 
predictions for the total cross section to be made, once I-induced topologies have been identified. 
Moreover, in order to reconstruct Q' and Wi, the separation of the current jet from the instanton 
band has to be as good as possible. We found the output of this procedure to be most useful in 
designing optimized observables in view of enhancing the separation power to normal DIS events. 
And finally, the reconstruction of Q' and £ would allow a boost to the instanton center of mass 
frame to be performed. 
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Figure 5: Reconstruction of Q' from the hadronic final state of I-induced events, a) Relative deviation 
of the true and the reconstructed value, b) Correlation between the mean value of the 



y2 



reconstructed Q and the true value. In both cases the upper figure is without and the lower 
figure is with a px cut on the current jet. 



The reconstruction of Q' is based on the assumption that the colour forces between the 
current jet and the rest of the partons, as modeled in the hadronisation phase, still allow the 
reconstruction of the current jet 4-vector. So q' = q — q" (cf. Fig. [[]) and the reconstruction of 
Q' is performed by reconstructing the four-vector q" of the particles belonging to the current 
jet. To identify the current jet the cone algorithm |2]| [26[] is applied in the hCMS. A cone of 
radius R = 0.5 (in the 77*-<J>-plane) turns out to perform best |H|. Motivated by the relative rf 
position of the current quark and the instanton band (cf. Fig. f|), we additionally require the jet 
to fulfill if > rj*, with rj* being the Et- weighted mean 77* value of all hadrons. The jet with the 
highest pt is then assumed to be the most likely candidate for the current jet. Fig. || shows that 
this procedure gives a rather good reconstruction of Q' , especially, if the jet is required to have 
a pt > 5 GeV. This is also a reasonable choice in order to have a clear signal in the detector. 
Once the current jet has been found, the hadrons belonging to the jet are removed from the final 
state and rj* is recalculated with the remaining hadrons in order to get a good estimator for the 
center of the instanton band. For the calculation of W, the four- vectors of all hadrons in a region 
left and right of 77* are considered. A band of 77* zb 1.1 gives the best possible results for the 
majority of events. This result is in agreement with the mean bandwidth visible in Fig. |] as well 
as with the expectation for an isotropic event. The reconstruction of Wi is not perfect, though. 
If the bandwidth is determined for every single event (based on the known value of Wi), a 
rather broad distribution is obtained, which peaks only weakly at ±1.1. In many events, slightly 
different bandwidths would be required for an optimal reconstruction of Wi. The procedure 
using the averaged bandwidth results in a quite large error of the reconstructed quantity of Wi 
of approx. 25 % (for Pr(jet) > 5 GeV). 

The variable x' can be calculated once Q' and Wi have been reconstructed. However, the errors of 
the Q' and the Wi reconstruction interact in such an unfavourable manner that a reconstruction 
of x' seems not very meaningful. Finally, £ can be reconstructed via s, which comprises the 
four- vectors of all the hadrons used for the reconstruction of Q' and Wi. The reconstruction of 
£ works relatively well in a broad range of the variable. However, a reconstructed pair of q' and 
£ leads to inconsistencies when trying to perform a boost to the instanton center of mass system 
in 50% of the cases, so that either variable would have to be corrected. 



5 Search Strategies 

A search strategy for I-induced processes is naturally based on the characteristic properties of 
the hadronic final state. The goal is to isolate - based on Monte Carlo predictions - an instanton 
enriched data sample by applying cuts on a set of observables. More than 30 observables have 
been investigated in |TI|[]. Fig. |S| demonstrates the procedure by which we measure the "sepa- 
ration power" of each observable. As an example, we investigate the distribution of the number 



6 For the results shown here additional cuts have been introduced to get a more realistic representation of the 
typical acceptance of a HERA detector. For this reason, some of the results presented here slightly deviate from 
those given in Q . 
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Figure 6: a) Shape normalized distributions of the number of charged particles in the instanton band 
n' B (without the hadrons belonging to the current jet) for normal DIS and I-induced events, 
b) Separation power (:= INS e ff/DIS e ff) after applying cuts in n' B for every possible cut. 
Explicitly shown are the values for the cut with the highest separation power, which fulfills 
the requirement of a remaining instanton efficiency of 10 %. 

of charged particles n' B , where the index "5" indicates that we consider only particles of the 
instanton band (77* ± 1.1), and the prime indicates that the hadrons belonging to the current jet 
are not considered. Starting from shape normalized distributions for normal DIS and I-induced 
events (Fig. ||a), we calculate the ratio of efficiencies for I-induced over normal DIS events for each 
possible cut value (Fig. []b). We require a minimum instanton efficiency, which we (arbitrarily) 
set to 10 %. The cut n' B > 12 leads to a ratio of the efficiencies, or - in our terminology - to 
a separation power of INS e ff/DIS e ff — 17. In general, applying cuts in only one observable 
typically leads to a separation power not higher than O(20) (c.f. [0). 

An obvious improvement in separating I-induced from normal DIS events is obtained by the 
combination of cuts in several observables. However, since most of the observables are highly 
correlated, the naive combination of the best one-dimensional cuts quite often fails in enhancing 
the separation power. On the contrary, in some cases it would even lead to a reduction of the 
separation power, an effect for which we give an example now. Since in any I-induced event at 
least one ss pair is produced, an excess in the number of neutral kaons in comparison to normal 
DIS events is expected. Fig. [7|a shows, however, that this excess is almost not visible for events 
with higher particle multiplicities. Investigating the correlation between the mean number of 
kaons and the number of charged particles (Fig. |7|b), we find the expected excess for relatively 
low particle multiplicities, only. For some value of n' B , the number of neutral kaons produced 
in normal DIS events actually exceeds the number produced in I-induced events. Applying cuts 
requiring a high value of n' B together with a large number of neutral kaons would therefore reduce 
the separation power. Since in addition the number of kaons produced rather strongly depends on 
variations of the underlying MC model, kaons play no role in our multi-dimensional cut-scenario. 
The study of strange mesons or baryons (notably A's) could, however, play a major role once an 
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instanton enriched sample is experimentally isolated from the data. 
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Figure 7: a) Shape normalized distributions of the number of neutral kaons in the instanton band for 
normal DIS and I-induced events, while applying the additional cut n' B > 7. b) Correlation 
between the mean number of neutral kaons and the number of charged particles. 

Based on a study of the correlations and mutual influences of different cuts, we now introduce 
a set of six observables which are able to enhance the separation power from O(20) to O(130). 
Fig. |8| shows these observables together with the cuts applied in each observable, indicated by the 
lines and the corresponding arrows^. 

This cut-scenario has been established in a study of what we call "reference" (or default) Monte 
Carlos. For the (default) simulation of normal DIS events we use the Lund MC generator 
ARIADNE ( ||27|| , including Pomeron exchange), which includes the Colour Dipol Model (CDM, 
p8| , [29| , ftl| ) to describe higher order perturbative QCD radiation. Although not perfect, 
ARIADNE presently gives the best description of the properties of the hadronic final state at 
HERA [[H], |32|. For the simulation of I-induced events we use QCDINS || |TB[, which is based 
on the matrix element of the I-subprocess and interfaced by default to HERWIG (T^j for the 
simulation of parton showers, fragmentation and particle decays (see section ||]). The shaded 
band shown in Fig. |8] stems from variations of the Monte Carlo simulations (see section ^)Q. 
All distributions presented here are normalized to the number of events. Fig. ||a shows the 
distribution of the pt of the current jet, i.e. the jet with the highest pt, for which rf(jet) > rf 
holds. The distribution for I-induced events is rather flat for p? values of 1 to ~ 7 GeV, while for 
normal DIS events the spectrum steeply falls towards higher p?. Fig. |8|b shows the distribution 
of Q' rec , a quantity which has no direct physical interpretation for normal DIS events, but 



which exhibits a rather good separation power (especially in correlation to n' B , c. f. \T2~\). The 



7 An explicit listing of these cuts is given in Figs. h3 and O. 

8 A list of all variations used in creating the shaded band in Fig. M is given in Fig. [12 and Fig. pi 
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Figure 8: Distributions of various observables for normal DIS and I-induced processes. Shown are 
the distributions for the "reference Monte Carlos" (INS markers = QCDINS + HERWIG, 
DIS markers = ARIADNE, including Pomeron exchange) and their variations (shaded band) 
resulting from the choice of different models or the variation of parameters of a model (c. f. 



Fig. 12 and Fig. 13). The lines and the corresponding arrows show the cut applied in each 
of the observables, with the arrows pointing in the direction of the allowed region. 



distribution of the transverse energy and of the charged particles in the instanton band (77* ±1.1) 
is displayed in Fig. |8]c and in Fig. [8]d. Both distributions reflect the high partonic (and following 
hadronic) activity in the I-subprocess. The observables depicted in Fig. ||e and [8]f make use 
of the shape of I-induced events which is supposed to be more isotropic than that of normal 
DIS events. Fig. ||e gives the distribution of the relative E in -E out difference, a quantity, which 
measures the Et weighted $ isotropy of an event (for a detailed description see |12|, |33|). The 
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more homogeneous the transverse energy is distributed in <3>, the lower the relative E in -E out 
difference becomes. Again, we investigate this quantity in the band[], after subtracting the 
hadrons belonging to the current jet. Normal DIS events seem to exhibit a surprisingly low value 
of this quantity, implying a rather isotropic $ distribution of the remaining hadrons. However, 
this behaviour strongly varies with the MC model used as is illustrated by the shaded band. 
This similarity of the distributions for normal DIS and I-induced events disappears for events 
with high Et as expected and as demonstrated in Fig. ^a. Here, normal DIS events exhibit the 
expected jet-like structure with low ($) isotropy, while I-induced events remain as isotropic as 
before. Finally, in Fig. B, we investigate the distributions of the first Fox- Wolfram moment 
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Figure 9: a) Shape normalized distributions of the relative Ei n -E ou t difference (in the instanton band, 
without jet hadrons) for normal DIS and I-induced events, while applying the additional cut 
E't b — 15 GeV. b) Correlation between the mean first Fox- Wolfram moment H\q and the 
number of charged particles n' B (in the instanton band, without jet hadrons). 

H\q ]34]] (normalized to the zeroth moment), an event shape variable which is independent of 
the axis chosen to calculate it. The Fox- Wolfram moments, in general, seem to perform better 
than all other studied event shape variables (like sphericity, thrust etc) in separating normal 
DIS from I-induced events (c.f. |12]). To understand, why the cut indicated by the arrow in 
Fig. [8]f improves the separation power, the correlation between H\q and n' B has to be taken into 
account (Fig. ||b). For high charged particle multiplicities, the mean value of H\q of normal 
DIS events lies 6e/ow;that of I-induced events so that cutting at higher values of H w is reasonable^. 

The multi-dimensional cut-scenario indicated in Fig. |8] leads to a separation power of 
INS e ff/DIS e ff ~ 130, with the remaining efficiencies INS eff ~ 10% and DIS eff ~ 0.08%. 

9 The separation obtained with this quantity can be enhanced by taking into account all hadrons in a broader 
range in rj* (rj* ± 1.6). This is indicated by the index "B*". 

10 The peak of the H\o distribution at H\q ~ 1 (Fig. |pf) is due to diffractive events (simulated by the Pomeron 
exchange or the SCI mechanism) which usually have very low multiplicities (n < 4). 
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If we extrapolate the prediction for the cross section given in Eq. ([/]) to the kinematical range 
investigated in this report (xbj > 10~ 4 ,Q 2 > 5 GeV 2 ,0.1 < y < 0.6), we find o HERA ~ 215 pb, 
and <J HERA — 22 pb after applying all cuts, respectively. For £ ~ 30 pb -1 , one finds for the 
number of I-induced events Nins — 670, while Ndjs — 1810 of normal DIS events are expected. 
This implies an (9(13) a effect when compared to the merely statistical DIS expectation. 



6 Dependencies on MC models 

In this section, we investigate how the final-state observables depend on the MC models used for 
the simulation of I-induced as well as normal DIS events. In the default version of QCDINS the 
cluster fragmentation model [|35| is used. Given the large number of O(10) partons produced in 



a narrow ?7*-interval by the I-subprocess, a thorough investigation of hadronization alternatives 
would be very desirable, since so far there is little experience with such configurations of high 
parton densities. In the following study we start with what is available in this context. For a 
selection of important distributions, Fig. pi] illustrates the effects of using the Lund string model 



22] , p6l p7jl as implemented in JETSET |2^| for fragmentation. The markers show the default 
implementation of QCDINS, while the shaded band reflects variations of QCDINS + JETSET, 
arising from different tunings of the JETSET parameters, which are obtained from studying the 



hadronic final state in e + e~ collisions ||38|| . Fig. |T0| a shows that the mean number of charged 
hadrons produced in an I-induced event slightly increases when using string fragmentation. For 
the number of neutral kaons (cf. Fig. [TD[b), however, a significant difference to the values obtained 
in the cluster fragmentation model is found. The Figs. |T0| c and d basically are a direct result of 
the larger number of hadrons produced, i.e. a slightly larger mean value of the transverse energy 
and a slightly more ($) isotropic distribution of the particles produced in an I-induced event. 
In general, we find that the prediction for the properties of the hadronic final state of I-induced 
events depends only weakly on the choice of the model used to simulate the hadronisation. This 



statement is supported by studying the effects of the HERWIG tuning, as suggested in plfl . Two 
parameters of the cluster fragmentation model influencing the maximum allowed cluster mass 
(CLMAX) and the mass distribution of the split clusters (PSPLT) are changed in the tuned 
HERWIG version. In Fig. |TI] we find, that this modifications strongly influence the distributions 
for normal DIS events^. On the contrary, the distributions of I-induced events only show rather 
slight variations^ Besides the influence of the p? of the current jet, the hadronic final state of 
I-induced events seems to be mainly determined by the dynamics of the I-subprocess, and by the 
available phase space in x' and Q' . 

Besides the parameters or models determining the hadronisation, we also varied the structure 
function of the proton from CTEQ4L |39| (which we used as default) to MRSH |3(J in the 



simulation of I-induced events. This variation results in a small decrease of the total cross section, 



n An exception is the number of neutral kaons, which seems to be not influenced by this parameter tuning. 

12 One exception to this general behaviour is the distribution of the first Fox- Wolfram moment, presented in 
Fig. gf, which shows rather strong variations when using JETSET or the untuned HERWIG instead of the tuned 
HERWIG for the hadronisation of I-induced events. 
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Figure 10: Comparison of different models to simulate the hadronisation of I-induced events with re- 
spect to their effect on a selection of distributions. The markers represent the default 
QCDINS implementation (QCDINS + HERWIG (tuned) = cluster fragmentation), while 
the shaded band represents QCDINS plus Lund string fragmentation (as implemented in 
JETSET) with different tunings applied (see text). 

but shows almost no visible effect in the distributions of the hadronic final state observables. 
Finally, we varied the scale ALL within its experimental uncertainty as explained in section |3|. 
Besides the expected strong influence on the cross section of I-induced events, the lower value of 

(3) 

A^- has the effect of slightly enhancing the charged particle multiplicity, while the higher value 
leads to slightly lower multiplicities. This effect can be easily explained as an increase or decrease 
in the mean gluon multiplicity which is proportional to l/a s and thus increases with decreasing 

A®. 

MS 

Fig. 



Yl shows the influence of all the variations of the simulation of I-induced events on the 
separation power of our multi-dimensional cut-scenario (explicitly given on top of the figure). For 
the simulation of normal DIS events we use (our default) MC model ARIADNE (with Pomeron 
exchange) as reference. As mentioned, the application of all cuts leads to a remaining DIS 
efficiency of 8 ■ 10 -4 or, respectively, a number of events of Njjjs — 1810, when normalized to 
C ~ 30 pb~ . The variations used in the modeling of I-induced events are sorted according 
to the separation power. The corresponding efficiencies range from approx. 9% to 11%. After 
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Figure 11: Comparison of the influence of HERWIG tuning for normal DIS and I-induced events with 
respect to its effect on a selection of distributions. 

normalization to the luminosity given above, one finds that the numbers of I-induced events 
remaining vary only slightly between Nins — 550 and Ni^s — 670. Because of its strong 
influence on the cross section for I-induced events, the obvious exception to this is the variation 
of the scale A™, which results in the worst and the best prediction for the remaining number of 
events, namely N INS — 140 and N INS ~ 1930, respectively. 

Another important aspect to investigate is the dependence of the separation power on variations 
of the MC models used to simulate normal DIS events. To this end, three different MC generators 



are used in our study: ARIADNE, LEPTO |4l[ and HERWIG. All generators are operated in 
the "tuned version" , i.e. with parameter settings obtained in optimizing the description of the 
hadronic final state at HERA [31|. For the generators ARIADNE and LEPTO, the simulation of 



diffractive events (modeled via the Pomeron exchange and the SCI-mechanism, respectively) has 
each been turned on and off. Fig. [H| shows the influence of these variations on the separation power 
of our multi-dimensional cut-scenario. This time, we fix the simulation of I-induced events to be 
described by our reference MC model QCDINS + HERWIG (tuned). The separation powers using 
LEPTO or ARIADNE turn out to be rather similar, while HERWIG gives a significantly different 
result. Quite remarkably, the expected difference in separation power, arising from the usage of 
different models for hadronisation, is even more pronounced in the tuned version of HERWIG 
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Figure 12: Dependence of the separation power (INS e ff/DIS e ff) of a multidimensional cut-scenario 
on the variation of MC models and parameters used to simulate I-induced events. The 
corresponding efficiencies are listed, as are the numbers of events remaining (assuming an 
integrated luminosity of £ ~ 30 pb~ ) . 



(which was tuned according to the same distributions as ARIADNE and LEPTO, c.f. ||31|| ). To 
demonstrate this effect, the results obtained using the untuned HERWIG implementation are 
additionally presented in Fig. [H| The obtained efficiencies result (for C ~ 30 pb _1 ) in numbers 
of remaining normal DIS events ranging from Npjs — 930, at DIS e ff ~ 4 • 10 -4 when using 
HERWIG (tuned), to N DIS ~ 2025, at DIS eff ~ 9 ■ 10" 4 , when using ARIADNE (without 
Pomeron exchange). For I-induced events we find the already quoted value of Nins — 669. 
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Figure 13: Dependence of the separation power (INS e tf/DIS e ff) of a multidimensional cut-scenario 
on the variation of the MC generator used to simulate normal DIS events. The correspond- 
ing efficiencies are listed, as are the numbers of events remaining (assuming an integrated 
luminosity of £ ~ 30 pb~ ) . 

7 Conclusion 

The experimental discovery of instanton-induced processes in DIS at HERA would be a novel, 
non-perturbative manifestation of QCD, and therefore be of basic significance. Based on the 
characteristic hadronic final state of I-induced events, we found a set of six observables that exhibit 
a good "separation-power" to normal DIS events. In our best cut-scenario, normal DIS events 
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are suppressed by a factor of 8 ■ 10~ 4 , while 10% of I-induced events survive. For C ~ 30 pb _1 
and the nominal value of the theoretical cross section, this results in a predicted number of 
Njns — 670, for I-induced events, while Nr>is — 1810 of normal DIS events are expected. Within 
the "bandwidth" of variations considered, the systematic uncertainties arising from the modeling 
of I-induced events are surprisingly small. In this case, the structure of the hadronic final state 
seems to be mainly determined by the dynamics of the I-subprocess and by the available phase 
space in x' and Q' . On the contrary, the simulation of normal DIS events is found to be much 
more sensitive to the variation of model parameters. Thus, a better understanding of the tails of 
distributions for normal DIS events turns out to be quite important. 
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